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Abstract 

Purpose Habitat loss is a significant cause of biodiversity 
loss, but while its importance is widely recognized, there is 
no generally accepted method on how to include impacts on 
biodiversity from land use and land use changes in cycle 
assessment (LCA), and existing methods are suffering from 
data gaps. This paper proposes a methodology for assessing 
the impact of land use on biodiversity using ecological 
structures as opposed to information on number of species. 
Methods Two forms of the model (global and local scales) 
were used to assess environmental quality, combining eco¬ 
system scarcity, vulnerability, and conditions for maintaining 
biodiversity. A case study for New Zealand kiwifruit pro¬ 
duction is presented. As part of the sensitivity analysis, 
model parameters (area and vulnerability) were altered and 
New Zealand datasets were also used. 

Results and discussion When the biodiversity assessment 
was implemented using a global dataset, the importance of 
productivity values was shown to depend on the area the 
results were normalized against. While the area parameter 
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played an important role in the results, the proposed alternative 
vulnerability scale had little influence on the final outcome. 
Conclusions Overall, the paper successfully implements a 
model to assess biodiversity impacts in LCA using easily 
accessible, free-of-charge data and software. Comparing the 
model using global vs. national datasets showed that there is 
a potential loss of regional significance when using the 
generalized model with the global dataset. However, as a 
guide to assessing biodiversity impact, the model allows for 
consistent comparison of product systems on an international 
basis. 

Keywords Biodiversity • Ecoregions • Impact assessment • 
Kiwifruit • Land use • Scale dependency 

1 Introduction 

In the drive for sustainability, knowledge of the environmen¬ 
tal impact of products and production systems is essential 
and life cycle assessment (LCA) is the prevailing framework 
for analyzing these impacts. Land use and land use changes, 
in particular extensive changes in agriculture or forestry, are 
also recognized as impacts that should be included in LCA 
but, despite being discussed since the early 1990s (Bare 
2010), no agreed method for this exists (Mila i Canals et al. 
2007). The relationship between land use and impacts on 
biodiversity has been long recognized (Barnes 1996) and 
loss of biodiversity is seen as one of the most severe threats 
to sustainability (Diaz and Cabido 2001). Despite recent 
focus due to the Convention on Biological Diversity (UNEP 
1992), there are few signs of reductions in the rate of loss of 
biodiversity (Butchart et al. 2010) and changes in land use 
are often the most important driver (Barnes 1996; Chapin 
et al. 2000; Sala et al. 2000; Henry et al. 2008; Haines-Young 
2009; Lenzen et al. 2009). However, it is not only major 
changes in land use that are important; small modifications to 
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land use that are not generally noted when land use changes 
are assessed, which might be equally important for biodiver¬ 
sity transformation (Barnes 1996; Haines-Young 2009). 
Barnes (1996) warns of the trap of putting higher value on 
things that exhibit immediate or obvious merit or which are 
easily measurable. He argued that a focus on rich ecological 
systems (e.g., tropical forests) is “inevitably a reactionary 
rather than proactive course of action” (p. 241) and that 
focusing on threatened hotspots can be ineffective. 

Assessment of land use using LCA is typically done as a 
product of the area affected, duration of the impact, and a 
qualitative measure of changes due to the impact (cf. 
Lindeijer 2000; Mila i Canals et al. 2007). Most proposals 
have used changes in species richness as a general measure 
for changes to biodiversity. The use of vascular plants as 
indicator has in particular been common in LCA (e.g., 
Muller-Wenk 1998; Kollner 2000; Weidema and Lindeijer 
2001; Toffoletto et al. 2007; Schmidt 2008; De Schryver 
et al. 2010), but other taxonomic groups are used as well 
(e.g., Geyer et al. 2010; de Baan et al. 2013). Use of selected 
taxonomic groups to represent biodiversity is debatable, 
partly because changes in species diversity only capture parts 
of the concept of biodiversity (cf. Millennium Ecosystem 
Assessment 2005) and partly because species richness in one 
taxonomic group has variable correlation with that of other 
taxonomic groups (e.g., Prendergast and Eversham 1997; 
van Jaarsveld et al. 1998; Lawton et al. 1998; Simila et al. 
2006). Some studies have shown correlations among taxo¬ 
nomic groups (e.g., Koellner 2002; Kati et al. 2004), but it is 
questionable whether these results are valid outside the area 
where the background datasets were collected. Interestingly, 
Gottfried et al. (2012) have documented that an early re¬ 
sponse to a long-term negative impact on biodiversity might 
in fact be a short-term increase in species richness, further 
stressing the fact that assessing changes in species richness 
can be an inadequate indicator for impacts on biodiversity 
since an increase cannot necessarily be interpreted as re¬ 
duced impact. Partly in response to the shortcomings of 
using species richness as a surrogate for biodiversity, some 
authors have proposed focusing on other features known to 
be important for biodiversity, such as water stress (Mila i 
Canals et al. 2009, 2010), structural features (Michelsen 
2008) or a framework based on expert opinions (Penman 
et al. 2010). Another problem with the proposed methods is 
the resolution. Some methods are specific to their case study 
areas and cannot be used elsewhere without intensive data 
gathering (e.g., Michelsen 2008; Geyer et al. 2010), some are 
assuming the same impact irrespective of where the land use 
occurs (e.g., Goedkoop et al. 2009), and yet others have a 
biogeographic subdivision, but on the scale of biomes (e.g., 
de Baan et al. 2013). Curran et al. (2011) pointed out that 
most approaches to biodiversity assessment in LCA are 
developed for specific regions, and criticized the global 


approach at biome level as being too coarse, but even at this 
scale, de Baan et al. (2013) experienced lack of data with a 
species approach. 

The main purpose of this paper is to present a model that 
enables assessments of land use impacts on biodiversity in 
life cycle impact assessment (LCIA), independent of any 
particular biogeographic region, based on globally available 
data with a higher resolution than biomes. The proposed 
method is exemplified with a case study for kiwifruit pro¬ 
duction in New Zealand. We also performed a sensitivity 
check for some of the most important assumptions in the 
model in order to test the strength of the method. 

2 Background 

Land cover can be classified at different scales. One classi¬ 
fication divides the world into 867 distinct ecoregions, which 
are defined as relatively large units of land containing dif¬ 
ferent natural communities and species, with boundaries that 
approximate the original extent of natural communities be¬ 
fore major land-use change (Olson et al. 2001). An interac¬ 
tive map of these ecoregions, as well as an indicator for the 
state of conservation for each, can be found on the World 
Wildlife Lund “WWL WildFinder” website (WWF 2012). 
Nine ecoregions can be found in New Zealand’s North and 
South islands (Fig. 1), all of them endemic to New Zealand. 
Ecoregions are the basis for our generic model, following 
Michelsen (2008). 

2.1 Land use in New Zealand 

New Zealand also has a national classification system, the 
Land Environments of New Zealand (LENZ). LENZ uses a 
combination of climatic and physiographic (landform and 
soil data) variables to classify New Zealand into areas con¬ 
taining ecosystems of similar type. LENZ is constructed 
hierarchically with four different resolutions, levels 1-4 di¬ 
viding New Zealand into 20, 100, 200, and 500 different 
environments, respectively (MfE 2009). The “Threatened 
Ecosystem Classification” system provides information on 
vulnerability, including the amount of remaining indigenous 
cover and protected areas for all 500 environments included 
in LENZ (Landcare Research 2012). We use level 4 of LENZ 
when we use our model specifically to assess kiwifruit 
production. 

In New Zealand, an estimated 50 % of land cover is 
indigenous vegetation, nearly 40 % is pasture, with the 
balance being exotic vegetation (9 %) and urban and built- 
up areas (0.8 %); only 1.5 % of land cover in New Zealand is 
horticulture (MfE 2007). Although large areas of New 
Zealand are protected for conservation purposes, Walker 
et al. (2005) claim that protection is skewed by a tendency 


<£) Springer 



Int J Life Cycle Assess (2014) 19:285-296 


287 


Fig. 1 Ecoregions of New 
Zealand and their state of 
conservation, based on 
information from the World 
Wildlife Fund's “WWF 
WildFinder” (WWF 2012) 
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to protect more wet mountainous environments while habi¬ 
tats and ecosystems in productive lowland and montane 
environments are less protected. As a consequence, indige¬ 
nous vegetation is mostly found in the alpine and upper 
montane zones and only traces on warmer, lower montane 
and lowland zones where productive land-use activities are 
more widespread. 

2.2 Kiwifruit production case study 

In this paper, we use kiwifruit production in New Zealand to 
illustrate a method for assessing impacts from land use on 
biodiversity. Our focus here is limited to impacts on biodi¬ 
versity caused by land use and we do not include impacts 
from land use changes. Because the land use impacts of an 
agricultural product’s life cycle are expected to be dominated 


by the cultivation phase, land use impacts from other life 
cycle stages are not included in this study. 

New Zealand is the second largest producer of kiwifruit in 
the world (Kilgour et al. 2008; FAO 2012). In 2010, the New 
Zealand total production was 378,5001, constituting 28 % of 
world production and a harvested area of ~ 12,800 ha (FAO 
2012). New Zealand’s kiwifruit industry has a high focus on 
environmental performance including efforts in carbon man¬ 
agement, use of water resources, lean manufacturing, pack¬ 
aging, transportation and carbon footprint assessment 
(Peltzer et al. 2011; ZESPRI 2010, 2012). 

In order to cover most of the areas used for kiwifruit 
production in New Zealand, 10 major kiwifruit growing 
regions are included in this study (Auckland, Gisborne, 
Hawke's Bay, Katikati, Nelson, Northland, Tauranga, Te 
Puke, Waikato, Whakatane). These orchards are situated in 
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only three ecoregions: Northland temperate kauri forests; 
North Island temperate forests; and Richmond temperate 
forests (Fig. 1). In total, kiwifruit production in the selected 
regions occupies 8,900 ha. 

As a starting point, we base our model on a method 
proposed by Michelsen (2008) for boreal forest, adapting it 
for more generic use by including aspects of a method 
proposed by Brentrup et al. (2002). For our kiwifruit case 
study, the functional unit is 1,000 tray equivalents (TE, 1 tray 
equivalent corresponds to 3.6 kg of green kiwifruit at gate), 
using data for the production year 2008/09 as an example. 


3 Methods 

3.1 Method description 

According to Mila i Canals et al. (2007), three aspects must 
be assessed when quantifying land use impacts on biodiver¬ 
sity: the area affected must be identified, the duration of the 
impact must be assessed, and a quantitative measure of 
impact on biodiversity must be established and assessed. 

Michelsen (2008) suggested an indirect approach for as¬ 
sessment of biodiversity quality ( Q ) and changes in quality 
at a given location based on ecosystem scarcity (ES), eco¬ 
system vulnerability (EV), and conditions for maintained 
biodiversity (CMB), using the formula 

g = ESxEVxCMB. (1) 

“Ecosystem scarcity,” originally proposed for LCIA by 
Weidema and Lindeijer (2001), is a measure of the intrinsic 
rareness of an ecosystem. In this paper, we will apply the ES 
equation proposed by Michelsen (2008); 

ES = I-A2L (2) 

A max 

where A pot is the potential distribution of the structure in 
focus and A max is the potential distribution of the most 
widespread structure at the relevant level. By the term “struc¬ 
ture,” we mean an ecologically and geographically defined 
area separated from other areas due to species composition 
and/or geophysical conditions. The scale-independent-term 
structure is used since different spatial scales can be used: 
biomes, ecoregions, etc. Michelsen (2008) suggested the use 
of ecoregions since these are the globally available data with 
the highest resolution (Olson et al. 2001; WWF 2012). Note 
that in this model, the structures with the highest scarcity 
thus get a value close to 1. 

“Ecosystem vulnerability” is a measure of the present 
condition of the structure. Different approaches are sug¬ 
gested (Peter et al. 1998; Weidema and Lindeijer 2001; 
Michelsen 2008). In WWF’s (2012) WildFinder, these are 


given on a three-grade ordinal scale. In order to translate this 
into a globally available quantitative scale, Michelsen (2008) 
suggested that the value 1.0 be used for ecoregions with a 
critical conservation status, 0.5 for ecoregions with a vulner¬ 
able conservation status, and 0.1 for intact ecoregions. 

“CMB” is an index for the actual impact on biodiversity in 
the affected area, ranging from 1, which represents no impact 
on biodiversity, to 0, which represents a complete removal of 
biodiversity. The CMB index is constructed from a suite of 
indicators known to be important for biodiversity in an 
ecosystem, and will thus give an indirect measure of the 
impact on biodiversity (cf. Larsson 2001). Michelsen 
(2008) constructed the CMB index for boreal forests using 
indicators such as amount of dead wood and alien species. 
Such indicators are not directly transferrable to other ecosys¬ 
tems and a challenge with this method is the need to con¬ 
struct ecosystem-specific indexes. To overcome this chal¬ 
lenge, we propose the use of hemeroby values (levels of 
naturalness) as suggested by Brentrup et al. (2002). 
Hemeroby values measure the human influence on ecosys¬ 
tems by determining the deviation from naturalness as a 
result of a specific land use (Kowarik 1999; Brentrup et al. 
2002; Walter and Stutzel 2009). The scale proposed by 
Brentrup et al. (2002) provides the naturalness degradation 
potential (NDP) for 39 land-use types. The scale includes 
most traditional uses of the land and covers most land uses 
across the world. 

Paracchini and Capitani (2011), in their summary of 
hemeroby, state that the degree of hemeroby is the result of 
the impact on a particular area and the organisms which 
inhabit it. The “cultivation of special crops (e.g., fruits vine)” 
or “extensive arable land use type” both have NDP of 0.7 (cf. 
Tables 1 and 2 from Brentrup et al. 2002). The 0.7 value is 
used as a proxy for kiwifruit orchards and when fitted to the 
scale of 0 for worst to 1 for best case; in this case study, the 
CMB is 0.3 for all orchards. 

The impact of land use is calculated as the change in 
quality as a result of the ongoing activity. Both Mila i Canals 
et al. (2007) and Michelsen (2008) consider a continued land 
use, such as an orchard, as a suppression of a relaxation 
process, i.e., as a delay in recovery from an impact to eco¬ 
system quality. In other words, the assumed alternative to the 
assessed land use is that the area recovers to its potential 
vegetation. The change in quality can then be calculated as 


AQ = Q t pot-fti, (3) 


where Q tpot is the potential quality of the area (defined as 
ES x EV) and Q n the quality of the area at the time (A) of the 
activity (calculated as ESxEVxCMB^). As a consequence, 
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Table 1 Land use impact on 
biodiversity from kiwifruit plan¬ 
tations following an adjusted 
version of the method proposed 
by Michelsen (2008) using 
hemeroby as CMB values 


TE kiwifruit tray equivalents, 

A pot potential distribution, ES 
ecosystem scarcity, EV ecosys¬ 
tem vulnerability, CMB condi¬ 
tions for maintained biodiversity, 
Q biodiversity quality 


Location 

TE/ha 

A pot 

ES 

A Q 

ha/y 

A(9 x haxy 

Ecoregion 

EV 

Nelson 

6,509 

13,200 

0.997 

0.349 

0.154 

0.054 

AA0408 

0.5 (vulnerable) 

Te Puke 

8,451 

84,400 

0.982 

0.687 

0.118 

0.081 

AA0405 

1 (critical) 

Tauranga 

8,125 

84,400 

0.982 

0.687 

0.123 

0.085 

AA0405 

1 (critical) 

Katikati 

7,411 

84,400 

0.982 

0.687 

0.135 

0.093 

AA0405 

1 (critical) 

Auckland 

7,409 

29,900 

0.994 

0.695 

0.135 

0.094 

AA0406 

1 (critical) 

Whakatane 

6,040 

84,400 

0.982 

0.687 

0.166 

0.114 

AA0405 

1 (critical) 

Northland 

5,839 

29,900 

0.994 

0.695 

0.171 

0.119 

AA0406 

1 (critical) 

Waikato 

5,770 

29,900 

0.994 

0.695 

0.173 

0.121 

AA0406 

1 (critical) 

Gisborne 

5,173 

84,400 

0.982 

0.687 

0.193 

0.133 

AA0405 

1 (critical) 

Hawke's Bay 

4,868 

84,400 

0.982 

0.687 

0.205 

0.141 

AA0405 

1 (critical) 


AQ = Q m -Qtx = (ES X EV)-(ES X EV x CMBq) 

= (1-CMB,!)(ES X EV). (4) 

3.2 Sensitivity of the method 

This section addresses the sensitivity of our method assump¬ 
tions, because (1) it is possible to have different perspectives 
when choosing A max ; (2) different spatial resolution can be 
selected for the structures in focus; and we discuss (3) the 
choice of values for EV in Michelsen (2008) and consider 
other values. 

3.2.1 Choice of A max 

Michelsen (2008) normalized the value of ES using the most 
widespread ecoregion as A max , therefore using the Sahara 


Desert which has a distribution of 4,639,900 km (WWF 
ecoregion PA 1327) as A max . Here, we test the use of the most 
widespread ecoregion in New Zealand, AA0405 North Is- 
land temperate forests with a distribution of 84,400 km , and 
the sum of the distribution of all ecoregions in New Zealand, 
in total 265,000 km , as alternatives. When A max =A pot , ES 
will necessarily be 0 and all impacts in such an ecoregion 
will by definition be 0. To avoid this, A max is set to 
90,000 km 2 when AA0405 is used as A max . 

3.2.2 Use of national classifications 

While ecoregions divide the ecosystems of the world on a 
global scale, many countries have, at least for country- 
specific purposes, classifications with a higher resolution. 
Here, we use the level 4 classification of LENZ as an alter¬ 
native to ecoregions. Regarding the vulnerability of the 
environments included in our model, New Zealand also has 


Table 2 Land use impact on biodiversity from kiwifruit plantations when A max is changed from using values for ecoregion PA1327 (Sahara Desert) 
to total area of ecoregions in New Zealand (TotNZ) and the most widespread ecoregion in New Zealand (90,000 km ) 


Location 

ES 



A Q 



AQx ha x y 




A 

^ max 

PA 1327 

A 

^ max 

TotNZ 

A 

max 

90,000 km 2 

A 

^ max 

PA 1327 

A 

^ max 

TotNZ 

A 

^ max 

90,000 km 2 

A 

^ max 

PA 1327 

A 

^ max 

TotNZ 

A 

max 

90,000 km 2 

Nelson 

0.997 

0.950 

0.853 

0.349 

0.333 

0.299 

0.054 

0.051 

0.046 

Te Puke 

0.982 

0.682 

0.062 

0.687 

0.239 

0.022 

0.081 

0.028 

0.003 

Tauranga 

0.982 

0.682 

0.062 

0.687 

0.239 

0.022 

0.085 

0.029 

0.003 

Katikati 

0.982 

0.682 

0.062 

0.687 

0.239 

0.022 

0.093 

0.032 

0.003 

Auckland 

0.994 

0.887 

0.668 

0.695 

0.311 

0.234 

0.094 

0.042 

0.032 

Whakatane 

0.982 

0.682 

0.062 

0.687 

0.239 

0.022 

0.114 

0.040 

0.004 

Northland 

0.994 

0.887 

0.668 

0.695 

0.311 

0.234 

0.119 

0.053 

0.040 

Waikato 

0.994 

0.887 

0.668 

0.695 

0.311 

0.234 

0.121 

0.054 

0.041 

Gisborne 

0.982 

0.682 

0.062 

0.687 

0.239 

0.022 

0.133 

0.046 

0.004 

Hawke's 

Bay 

0.982 

0.682 

0.062 

0.687 

0.239 

0.022 

0.141 

0.049 

0.004 


ES is ecosystem scarcity, Q is biodiversity quality 
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the Threatened Ecosystem Classification, which is based on 
the assessment of indigenous cover. To get more accurate data, 
the Land Cover Database, LCDB2 (MfE 2009) was used in 
combination with LENZ level 4 to determine in which LENZ 
environments the orchards occur and how much of the indig¬ 
enous vegetation is left in those environments. QuantumGIS 
(QGIS 2012), an open-source geospatial information system 
software, was used for data processing. 

To identify where orchards are located, we used the most 
up-to-date, publicly available version of land cover data (the 
second edition of New Zealand’s Land Cover Database 
(LCBD2)) from the Koordinates portal (www.koordinates, 
com). The polygons of the LCDB2 class “orchards, 
vineyard, and perennial crops” were used as boundaries for 
the areas of each LENZ level 4 environment, and the area of 
each LENZ level 4 environment in the orchards was 
calculated. This process was carried out for each of the 10 
selected kiwifruit growing regions. Following this approach, 
the ecosystem vulnerability cannot be classified following 
the ordinal scale Michelsen (2008) originally used. Here, 
Eq. (6) (below) is used to calculate values for ecosystem 
vulnerability for each LENZ structure based on amount of 
indigenous cover left. 

LENZ has much higher resolution than the data on the 
kiwifruit cultivation areas, and more precise geospatial data 
for the location of orchards were unavailable. Therefore, ES 
was calculated using Eq. (5), which is a weighted average for 
each kiwifruit growing region, following 

K +((+)+ (5) 

where ES is the weighted average of the scarcity of the 
area in focus (as a kiwifruit growing region), A vot . is the 
potential distribution of environment i, A max is the most 
widespread environment at the same level (here LENZ 
level 4), and LU; is the percentage of environment i for 
the total area included. 

3.2.3 Scaling ecosystem vulnerability 

Michelsen (2008) used the state of conservation from 
WWF’s Wildfmder (WWF 2012) for scaling ecosystems' 
vulnerability. Here, this is adjusted to an exponential func¬ 
tion based on knowledge of the fraction of the original 
vegetation lost (FL). Following the idea from Michelsen 
(2008), EV=0 when the entire structure is still present 
(FL=0), and EV=1 when the entire structure is gone 
(FL=1). In order to construct the function, it was also as¬ 
sumed that when 50 % of the fraction was lost, EV is 0.25, 
giving the function 

EV =-0.125 +0.125exp(2.1972 x FL). (6) 


The steepness of the curve (Fig. 2) is debatable, and a 
universal curve probably does not exist (cf. de Schryver et al. 
2010). Lindgaard and Henriksen (2011) provide estimates on 
the relationship between fraction of structures lost and how 
endangered they are. Based on this information combined 
with Eq. (6), we suggest that 1.0 is still used for critical 
ecoregions, but that vulnerable ecoregions are adjusted to 
0.45 and intact ecoregions to 0.09. The values based on this 
approach are thus comparable with the values applied by 
Michelsen (2008). 

As for ecosystem scarcity, ecosystem vulnerability (EV) 
must also be assessed as a weighted average when more than 
one structure is included in the area considered, following 

n 

EV = ^2 (EV, X LU,), (7) 

/= 1 


where EV is the weighted average, EV; is the ecosystem 
vulnerability of structure i following Eq. (6), and LU; is the 
percentage of the total land use found in structure i. 


4 Results 

Table 1 shows the results of the land use impact of producing 
1,000 TE at the included locations following the proposed 
default values for ES, EV, and CMB, where hemeroby values 
are used for CMB. The results are sorted from lowest to 
highest impact. 

Figure 3 is a graphical representation of parameters of the 
model and of the results. Note that the values are adjusted to 
fit a 0 to 1 scale for graphical representation. 

Table 2 shows the results when A max is changed to the 
sum of all ecoregions in New Zealand (TotNZ) and to the 



0 0.2 0.4 0.6 0.8 1 


Ecosystem vulnerability 

Fig. 2 Graphical presentation of the relationship between fraction lost 
of a structure and its ecosystem vulnerability, following the proposed 
equation (Eq. (6)) 
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Fig. 3 Graphical presentation of 
the model parameters and results 
shown in Table 1 . * for TE/ha, 
the scores are in 1CT 5 . ** forT pot , 
the scores are in 10 6 



largest New Zealand area (here 90,000 km is used to avoid 
A max being equal to A pot ). 

The nonlinear adjustment of EV (see Section 3.2.3) has only 
minor implications when used at ecoregion level and for New 
Zealand kiwifiuit production is only relevant for the Nelson 
region, where EV is reduced from 0.50 to 0.45. When the score 
for Nelson shown in Table 1 is recalculated with this adjustment, 
the total score is reduced from 0.054 to 0.048. The choice of 
A max has a larger impact (Table 2); the absolute values of ES and 
A Q decline with a smaller A max (Table 2). This influences how 
the different aspects impact on calculating the total impact from 
land use. As can be seen from Fig. 4, when A max is reduced 
from the size of the Sahara (4,639,900 km ) to the area of New 
Zealand (265,000 km ) or the most widespread ecoregion in 


New Zealand (90,000 km ), not only are the absolute scores 
reduced, but the ranking among the kiwifruit growing regions 
also changes. 

In Table 3, the results using LENZ level 4 are shown; the 
underlying data showing the impact from different structures 
are given in the Electronic Supplementary Material. 

5 Discussion 

At present, there is no generally accepted method for how to 
include impacts from land use and land use changes in LCA. 
For land use impacts on biodiversity, several of the methods 
proposed make use of changes in species diversity and in 


Fig. 4 Total land use impact on 
biodiversity for the different 
locations whenT max is altered 
from the most widespread 
ecoregion worldwide (Sahara 
Desert), to the total area of 
ecoregions in New Zealand or 
the most widespread ecoregion 
in New Zealand. For Nelson*, 
EV is reduced from 0.50 to 0.45 
(see text) 
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Table 3 Land use impact on biodiversity from kiwifruit plantations 
using national data for assessing ES and EVon a detailed level. See text 
for more details 


Fo cation 

TE/ha 

ES 

EV 

CMB 

A Q 

ha/y 

A(9 x haxy 

Nelson 

6,509 

0.969 

0.888 

0.3 

0.603 

0.154 

0.093 

Te Puke 

8,451 

0.729 

0.661 

0.3 

0.337 

0.118 

0.04 

Tauranga 

8,125 

0.799 

0.746 

0.3 

0.417 

0.123 

0.051 

Katikati 

7,411 

0.823 

0.734 

0.3 

0.423 

0.135 

0.057 

Auckland 

7,409 

0.831 

0.799 

0.3 

0.465 

0.135 

0.063 

Whakatane 

6,040 

0.912 

0.917 

0.3 

0.586 

0.166 

0.097 

Northland 

5,839 

0.849 

0.623 

0.3 

0.370 

0.171 

0.063 

Waikato 

5,770 

0.715 

0.884 

0.3 

0.442 

0.173 

0.077 

Gisborne 

5,173 

0.981 

0.946 

0.3 

0.650 

0.193 

0.126 

Hawke's Bay 4,868 

0.974 

0.974 

0.3 

0.664 

0.205 

0.136 


TE kiwifruit tray equivalents, ES ecosystem scarcity, EV ecosystem 
vulnerability, CMB conditions for maintained biodiversity Q biodiver¬ 
sity quality 


particular vascular plant diversity (e.g., Muller-Wenk 1998; 
Kollner 2000; Toffoletto et al. 2007; Schmidt 2008; 
Goedkoop et al. 2009). Species diversity only captures part 
of the concept of biodiversity, and the correlation between 
vascular plant diversity and species diversity for other taxo¬ 
nomic groups is in most cases very low. When impacts on 
biodiversity from land use and land use changes are to be 
incorporated in LCA, there is a need for indicators that are 
ecologically sound and easy to measure. Because the rela¬ 
tionship between vascular plant diversity and biodiversity is 
arguable, the use of vascular plants as a surrogate for biodi¬ 
versity worldwide seems to only fulfill the need for an 
indicator that is easy to measure and even then still depends 
on actual knowledge of vascular species diversity, which is 
not available for all ecosystems in the world (de Baan et al. 
2013). 

Our intention has been to demonstrate an alternative ap¬ 
proach where impacts on biodiversity can be assessed with¬ 
out the need to use any taxa as indicators. The method this 
study is based upon was originally developed for boreal 
forests but is here adapted to be used in any other region of 
the world, and an example for kiwifruit orchards in New 
Zealand is presented. Michelsen (2008) introduced an 
ecosystem-specific index (CMB) and used this to assess the 
impact on biodiversity in the ecosystem in focus given the 
identified use of the area. Such CMB indexes must conse¬ 
quently be developed for all ecosystems and scaled for all 
relevant activities. This requires data that are often not read¬ 
ily available, which is also the case for kiwifruit production 
in New Zealand. This lack of data is overcome by using the 
generic values for ecosystem naturalness provided by 
Brentrup et al. (2002). With this adjustment, the methodolo¬ 
gy proposed by Michelsen (2008) becomes applicable for 


almost any land use worldwide, since Brentrup et al. (2002) 
provide naturalness values for almost all land use situations 
and Michelsen’s (2008) proposal makes use of globally 
available data on ecosystem scarcity and ecosystem vulner¬ 
ability at ecoregion level (WWF 2012). The results for the 10 
different kiwifruit growing regions of our case study are 
shown in Table 1. As long as productivity and location of 
an ecoregion are known, calculations as shown in Table 1 are 
straightforward. 

It is noteworthy that a major argument for introducing the 
CMB index is that it makes it possible to evaluate different 
management regimes differently since not just changes in land 
use, but also modifications and changes in management could 
be of major importance for biodiversity (Bames 1996; Haines- 
Young 2009). For example, Michelsen (2008) showed how 
even small changes to management regimes could change the 
results significantly for timber production in a boreal forest. 
By using generic data, as in this paper, this is not possible 
since all similar land use activities (such as kiwifmit produc¬ 
tion) in an ecoregion will result in equal impact, irrespective of 
management differences, with the only difference being de¬ 
pendent on productivity. We thus face a trade-off between the 
ability to apply the method globally without a need for new 
data and the ability to differentiate between different levels of 
impact between different management regimes. 

The sensitivity of central assumptions made by Michelsen 
(2008) is tested here. To calculate ecosystem scarcity using 
A max , Michelsen (2008) proposed using the most widely 
distributed ecoregion (the Sahara) as A max . This follows the 
same rationale as originally proposed by Weidema and 
Lindeijer (2001). However, the final score, and consequently 
the importance of ES for the final value for land use impacts, 
depends on the choice of A max . When lower values are 
chosen for A max , ES and A Q will also be reduced. This is 
shown in Table 2 and Fig. 4. This will again influence how 
the different aspects affect the calculation of total land use 
impact—as can be seen from Fig. 4, when A max is reduced 
from the size of the Sahara (4,639,900 km ) to the area of 
New Zealand (265,000 km ) or the most widespread 
ecoregion in New Zealand (90,000 km ). When using differ¬ 
ent values ofH max , not only are the absolute scores reduced, 
but the mutual ranking among the plantation areas is also 
changed. 

When the Sahara Desert is used as A max , EV plays an 
important role in the kiwifruit case study presented here and 
clearly separates the Nelson orchard from the others. What 
really differentiates the other orchards is the land require¬ 
ment, i.e., the kiwifruit yield (Table 1). With reduced values 
of A max? land requirements become less important and the 
ecoregion in which the orchards are located increases in 
importance. 

From our point of view, there is a need to balance the total 
area requirement to the potential quality of these areas, but 
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with present knowledge, it is not possible to decide upon 
what is the correct level. This is, at least partly, a normative 
question; i.e., what is most important—to reduce the area 
requirement or to protect specific areas? This simply is a 
weighting problem: the less weight that is given to quality 
aspects of the land used, the more weight will be put on area 
demand, which is directly correlated to yield, irrespective of 
how quality is assessed. 

At this stage, at least for international comparisons, our 
recommendation is to use ecoregions, using the Sahara as 
default A max . There are two main arguments for this. First, 
the problem of lack of data is eliminated since these data are 
available for a global scale. Second, this also ensures con¬ 
sistency, which at present is hard to accommodate when 
using and comparing different national datasets. In the kiwi- 
fruit case study presented here, the use of Sahara as A max 
gave small differences in ES for the relevant ecoregions in 
New Zealand; however, since all ecoregions in New Zealand 
are relatively small and endemic to the country, and two of 
them, hosting 9 of the 10 orchards, are regarded to be in a 
critical state of biodiversity loss, this makes sense from a 
biodiversity-conserving perspective. However, it is impor¬ 
tant to have in mind that size of ecoregions is not fixed once 
and for all; with new knowledge, new subdivisions might 
occur and also global systems at a finer scale might be 
available. 

Equally important is the choice of structural level. Both 
Michelsen (2008) and our study (Tables 1 and 2) focused on 
ecoregions. Again, obvious reasons for this are globally 
available data and the ability to make international compar¬ 
isons. The problem is that data, at least for EV, are rather 
crude at this level and a very large number of areas are 
termed critical and get the highest score. As a consequence, 
large areas are treated equally and this is most likely not 
correct under all circumstances. For example, even though 
an ecoregion can be regarded as being in a critical condition, 
there are still structures at lower levels that are not endan¬ 
gered (Lindgaard and Henriksen 2011). 

We used the LENZ classification as an alternative to 
ecoregions to provide a more precise picture of the actual 
status of the areas in focus. Here, we chose to use LENZ 
level 4, increasing the number of potential area units in New 
Zealand from nine ecoregions to 500 LENZ types. For most 
of the orchards presented here, the use of LENZ indicated a 
lower impact of kiwifruit production on biodiversity since 
EV was assessed to be lower in all but one case. For the Te 
Puke and Northland orchards, the reduction in EV was close 
to 40 %. Interestingly, EV increased for the Nelson region 
when a value more specific to that environment was used, 
further stressing limitations of EV at the ecoregion level. 

There were 75 LENZ structures relevant for our case 
study, and for example, the orchards around Napier occurred 
in 17 level 4 structures (see the Electronic Supplementary 


Material). However, when only 3 to 5 of the LENZ environ¬ 
ments, covering most of each kiwifruit region, were includ¬ 
ed, the results did not change much (results not shown). It 
might thus be possible to obtain similar results using consid¬ 
erable less data; however, we do not attempt to state the 
minimum number of structures to be included as this would 
probably be case sensitive. Our assessment was done using 
publicly available land cover data (LCDB2) from which we 
determined the amount of orchards in each region, and in 
which environment from LENZ level 4 they were located. 
The results are subject to (1) the data quality of the geospatial 
information used, such as temporal correlation between 
LCDB2 (which is based on satellite images dated from 
2002) and the reference year used in this research and (2) 
the assumption that the distribution of kiwifruit orchards is 
proportional to the distribution of orchard and perennial 
crops (from LCBD2) in that environment. It is thus assumed 
that the geographic and temporal correlations sufficiently 
fulfill the aim of the study, providing a general picture of 
areas that have been mapped as orchards and perennial crops, 
which are or could be kiwifruit plantations. 

Still, we do not recommend this as a general approach. It 
is important to bear in mind that LCA should be a tool 
supplementary to other tools. The aim should not be to make 
LCA so detailed that users could regard it as an alternative to, 
for example, site-specific assessments for area and species 
protection or detailed planning (e.g., of precise locations for 
new plantations). These questions are probably better dealt 
with using other planning tools. We have shown here that the 
present EV values at ecoregion level are probably too coarse. 
The only region not given the highest score in Table 1 
(Nelson) has an average value when local data are used 
(Table 3). The level of detail needed (ecoregion or LENZ 
level 4) for impact assessment of land use is still an open 
question. 

We also propose a new equation for calculating EV 
(Eq. (6)). Although this did not have much influence when 
results were assessed at ecoregion level (Fig. 4), the ap¬ 
proach was also used for calculating EV at LENZ level 4 
(Eqs. (6) and (7)) and should give a result that is more sound 
ecologically, but as mentioned in Section 3.2.3, the steepness 
of this curve can still be debated. A universal curve probably 
does not exist and a steeper curve would make the differ¬ 
ences between EV more pronounced, giving relatively lower 
values for structures with higher remaining indigenous 
cover. 

A drawback of the use of hemeroby values as a proxy 
for the CMB index, introduced by Michelsen (2008), is 
the loss of ability to delineate different management re¬ 
gimes for similar activities in similar areas. It is possible 
to foresee situations where farmers would want to in¬ 
crease (local) biodiversity within their orchard (e.g., have 
protected areas set aside). However, this would increase 
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the area requirement, so the total impact per functional 
unit might still not change significantly. In this paper, we 
acknowledge that there is an obvious trade-off between 
data availability and level of precision of the assessments, 
keeping the focus on the opportunity to apply the method 
globally. Another drawback of the hemeroby is the weak 
scientific basis for the values. The values are assumed to 
be “an integrative measure of the impacts of all human 
interventions on ecosystems” (Paracchini and Capitani 
2011), but empirical data justifying the scores are mostly 
lacking. Despite this, they are assumed to be sufficiently 
robust to be recommended by, for example, the European 
Commission (cf. Paracchini and Capitani 2011). Also, this 
is a trade-off. De Baan et al. (2013) have tried to make 
impacts factors for eight different land use classes on 
biome level (14 worldwide) and shown data deficits even 
at this level. Despite their weaknesses, hemeroby values 
do offer a framework for global assessments applicable to 
a much finer spatial scale. 

Another drawback of the method is the interpretation of 
the scores in ecological terms. The value is given as quality 
changes multiplied by area (in hectares) and time (in years). 
This means that the value 1 represents a total removal of 
biodiversity in the most valuable structure (ES and EV both 
1) over 1 ha for 1 year. The values given are relative to this 
and changes irrespective of changes in area, time, intrinsic 
quality and/or impact on biodiversity in the area in focus. As 
mentioned earlier, the scores will be influenced by the 
choices and assumptions made; here, exemplified with dif¬ 
ferent scores on ES as a consequence of choice of A max and 
different scores on EV as a consequence of knowledge of 
fraction lost of a structure. 

6 Conclusions 

The goal of this study was to propose a globally applicable 
method that enables assessments of impacts on biodiversity 
from land use, using kiwifruit production in New Zealand as 
a case study. This was done by combining two previously 
proposed methods, where naturalness data from Brentrup 
et al. (2002) were used as default values for the CMB 
index (conditions for maintained biodiversity) proposed by 
Michelsen (2008) and combined with data on ecosystem 
scarcity and ecosystem vulnerability. Scarcity and vulnera¬ 
bility are parameters that are intrinsic to the area under study 
while the CMB index relates to the land use intervention. 
The CMB index brings to the model the principle of the 
change of quality caused by the land use. The use of the 
hemeroby concept as a proxy value for CMB is indeed a 
simplification, but allows the model to be applied globally as 
long as the area requirement and location in an ecoregion are 


known (as demonstrated by the kiwifruit case study). If 
hemeroby's linear approach can limit the use of the model 
for being generic on the one hand, the concept is easy to 
apply and does not require any extra data other than yield and 
location on the other hand. 

The sensitivity of the method was assessed by varying 
dmax and by different approaches to calculating EV. This 
showed that present values for EV at ecoregion level are 
probably too crude to give a realistic picture of the onsite 
situation in specific cases. We believe our assessment carried 
out using LENZ data with a higher spatial resolution gives a 
more realistic representation of the actual vulnerability of the 
areas in focus, but agree that such data are not readily 
available for all locations. Different choices of A max give 
different weights to the different aspects of land use. 

Acknowledgments The article is based on research initially conducted 
by the author Carla Coelho while at Landcare Research and was support¬ 
ed by Government capability funding to Landcare Research from the 
Foundation for Research, Science and Technology. At present, Carla 
works as Geospatial Analyst in the Auckland Council, and the views 
expressed in this article are independent from her current employer. Ottar 
Michelsen had financial support from the Norwegian Bioenergy Innova¬ 
tion Centre (CenBio). The authors acknowledge interesting discussions 
with colleague Pascale Michel, input from William Lee, who recom¬ 
mended the use of LENZ for better adaptation to the New Zealand 
context, and Nancy Golubiewski. Christine Bezar and Bob Frame 
reviewed an earlier draft. The authors also acknowledge thoughtful 
comments from two anonymous reviewers. 


References 

Bare JC (2010) Life cycle impact assessment research developments 
and needs. Clean Technol Environ Policy 12:341-351 

Bames N (1996) Conflicts over biodiversity. In: Sloep PB, Blowers A 
(eds) Environmental policy in an international context. Environ¬ 
mental problems as conflicts of interest. Arnold, London, pp 217- 
254 

Brentrup F, Kiisters J, Lammel J, Kuhlmann H (2002) Life cycle impact 
assessment of land use based on the hemeroby concept. Int J Life 
Cycle Assess 7:339-348 

Butchart SHM, Walpole M, Collen B, van Strien A, Scharlemann JPW, 
Almond REA et al (2010) Global biodiversity: indicators of recent 
declines. Science 328:1164-1168 

Chapin FS III, Zavaleta ES, Eviner VT, Naylor RT, Vitousek PM, 
Reynolds HL et al (2000) Consequences of changing biodiversity. 
Nature 405:234-242 

Curran M, de Baan L, de Schryver A, van Zelm R, Hellweg S, Koellner 
T, Sonnemann G, Huijbregts M (2011) Toward meaningful end 
points of biodiversity in life cycle assessment. Environ Sci 
Technol 45:70-79 

De Baan L, Alkemade R, Koellner T (2013) Land use impacts on 
biodiversity in LCA: a global approach. Int J Life Cycle Assess 
18(6): 1216—1230 

De Schryver AM, Goedkoop MJ, Leuven RSEW, Huijbregts MAJ 
(2010) Uncertainties in the application of the species area relation¬ 
ship for characterisation factors of land occupation in life cycle 
assessment. Int J Life Cycle Assess 15:682-691 


<£) Springer 



Int J Life Cycle Assess (2014) 19:285-296 


295 


Diaz S, Cabido M (2001) Vive la difference: plant functional 
diversity matters to ecosystem processes. Trends Ecol Evol 
16:646-655 

FAO (2012) FAOSTAT. http://faostat3.fao.org/home. Accessed 2 
August 2012 

Geyer R, Stoms DM, Lindner JP, Davis FW, Wittstock B (2010) 
Coupling GIS and LCA for biodiversity assessments of land use. 
Int J Life Cycle Assess 15(5):454^467 
Goedkoop M, Heijungs R, Huijbregts M, De Schryver A, Struijs J, van 
Zelm R (2009) ReCiPe 2008. A life cycle impact assessment 
methods which comprises harmonised category indicators at the 
midpoint and the endpoint level. First edition. Report I: Charac¬ 
terisation. Available at http://www.leidenuniv.nl/cml/ssp/publications/ 
recipe_characterisation.pdf Accessed 13 April 2012 
Gottfried M, Pauli H, Futschik A, Akhalkatsi M, Barancok P, Alonso 
JLB et al (2012) Continent-wide response of mountain vegetation 
to climate change. Nat Clim Change 2:111-115 
Haines-Young R (2009) Land use and biodiversity relations. Land Use 
Policy 26S: 178-186 

Henry PY, Lengyel S, Nowicki P, Julliard R, Clobert J, Celik T et al 
(2008) Integrating ongoing biodiversity monitoring: Potential ben¬ 
efits and methods. Biodivers Conserv 17:3357-3382 
Kati V, Devillers P, Dufrene M, Legakis A, Vokou D, Lebrun P (2004) 
Testing the value of six taxonomic groups as biodiversity indica¬ 
tors at a local scale. Conserv Biol 18:667-675 
Kilgour M, Saunders C, Scrimgeour F, Zellman E (2008) The key 
elements of success and failure in the NZ kiwifruit industry. 
Research Report No. 311. Agrobusiness and Economics Research 
Unit, Lincoln, New Zealand 

Koellner T (2002) Land use in product life cycles and its consequences 
for ecosystem quality. Int J Life Cycle Assess 7:130 
Kollner T (2000) Species-pool effect potentials (SPEP) as a yardstick to 
evaluate land-use impacts on biodiversity. J Clean Product 8:293- 
311 

Kowarik I (1999) Naturlichkeit, Naturnahe und Hemerobie als 
Bewertungskriterien. In: Konold W, Bocker R, Hampicke U 
(eds) Handbuch Naturschutz und Landschaftspflege. Ecomed, 
Landsberg 

Landcare Research (2012) Threatened Environment Classification. 
http://www.landcareresearch.co.nz/resources/maps-satellites/ 
threatened-environment-classification, accessed 13 August 2012 
Larsson TB (ed) (2001) Biodiversity evaluation tools for European 
forests. Ecological Bulletins 50. Blackwell Science, Oxford 
Lawton JH, Bignell DE, Bolton B, Bloemers GF, Eggleton P, 
Hammond PM et al (1998) Biodiversity inventories, indicator 
taxa and effects of habitat modification in tropical forest. Nature 
391:72-76S 

Lenzen M, Lane A, Widmer-Cooper A, Williams M (2009) Effects of 
land use on threatened species. Conserv Biol 23:294-306 
Lindeijer E (2000) Biodiversity and life support impacts of land use in 
LCA. J Cleaner Prod 8:313-319 

Lindgaard A, Henriksen S (eds) (2011) Norsk rodliste for naturtyper 
2010. Artsdatabanken, Trondheim 

MfE, Ministry for the Environment (2007) Environmental reporting, 
land cover, Available at http://www.mfe.govt.nz/environmental- 
reporting/land/cover/. Accessed 25 November 2011 
MfE, Ministry for the Environment (2009) Publications, state of the 
environment, Land Environments of New Zealand, Available at 
http://www.mfe.govt.nz/publications/ser/lenz-apr03.html. 
Accessed 13 August 2012 

Michelsen O (2008) Assessment of land use impact on biodiversity. 
Proposal of a new methodology exemplified with forestry opera¬ 
tions in Norway. Int J Life Cycle Assess 13:22-31 
Mila i Canals L, Bauer C, Depestele J, Dubreuil A, Freiermuth 
Knuchel R, Gaillard G et al (2007) Key elements in a framework 


for land use impact assessment within LCA. Int J Life Cycle 
Assess 12:5-15 

Mila i Canals L, Chenoweth J, Chapagin A, Orr S, Anton A, Clift R 
(2009) Assessing freshwater use impacts in LCA: Part I - inventoiy 
modelling and characterisation factors for the main impact path¬ 
ways. Int J Life Cycle Assess 14:28^42 
Mila i Canals L, Chapagin A, Orr S, Chenoweth J, Anton A, Clift R 
(2010) Assessing freshwater use impacts in LCA, part 2: case 
study of broccoli production in the UK and Spain. Int J Life Cycle 
Assess 15:598-607 

Millennium Ecosystem Assessment (2005) Ecosystems and human 
well-being: current state and trends: findings of the Condition 
and Trends Working Group / co-chairs Hassan R and Scholes R. 
Island, Washington, DC 

Muller-Wenk R (1998) Land use—the main threat to species. How to 
include land use in LCA. IWO—Diskussionsbeitrag 64. Institute 
for Economy and the Environment (IWO). University St. Gallen, 
St. Gallen 

Olson DM, Dinerstein E, Wikramanayake ED, Burgess ND, Powell 
GVN, Underwood EC et al (2001) Terrestrial ecoregions of the 
world: a new map of life on Earth. BioScience 51:933-938 
Paracchini ML, Capitani C (2011) Implementation of a EU wide indi¬ 
cator for the rural-agrarian landscape. JRC Scientific and Techni¬ 
cal Reports EUR25114EN-2011. European Commission, Joint 
Research Centre, Brussels 

Peltzer DA, MacLeod CJ, Gormley AM, Perry M, Burrows L, Moller H, 
Benge J (2011) Weed risk, plant biodiversity and carbon storage 
within kiwifruit orchard land titles. Landcare Research Contract Re¬ 
port LC626 to ZESPRI International Ltd and Bay of Plenty Regional 
Council, New Zealand. Landcare Research, Dunedin, New Zealand 
Penman T, Law B, Ximenes F (2010) A proposal for accounting for 
biodiversity in life cycle assessment. Biodivers Conserv 19:3245— 
3254 

Peter D, Krokowski K, Bresky J, Petterssn B, Bradley M, Woodtli H, 
Nehm F (1998) LCA graphic paper and print products (part 1, long 
version). Infras AG (Zurich), Axel Springer Verlag AG 
(Hamburg), Stora (Falun, Viersen) and Canfor (Vancouver) 
Prendergast JR, Eversham BC (1997) Species richness covariance in 
higher taxa: empirical tests of the biodiversity indicator concept. 
Ecography 20:210-216 

QGIS, Quantum GIS project, version 1.7.4 Wroclaw (2012) http:// 
www.qgis.org. Accessed 5 March 2013 
Sala OE, Chapin FS III, Armesto JJ, Berlow E, Bloomfield J, Dirzo R 
et al (2000) Global biodiversity scenarios for the Year 2100. 
Science 287:1770-1774 

Schmidt JH (2008) Development of LCIA characterisation factors for 
land use impacts on biodiversity. J Cleaner Product 16:1929-1942 
Simila M, Kouki J, Monkkanen M, Sippola AL, Huhta E (2006) Co¬ 
variation and indicators of species diversity: can richness of forest¬ 
dwelling species be predicted in boreal forests? Ecol Indicat 
6:686-700 

Toffoletto L, Bulle C, Godin J, Reid C, Deschenes L (2007) LUCAS—a 
new LCIA method used for a Canadian-specific context. Int J Life 
Cycle Assess 12:93-102 

UNEP, United Nations Environment Programme (1992) Convention on 
Biological Diversity. Text and annexes. UNEP, Geneva 
van Jaarsveld AS, Freitag S, Chown SL, Muller C, Koch S, Hull H et al 
(1998) Biodiversity assessment and conservation strategies. Sci¬ 
ence 279:2106-2108 

Walker S, Price R, Rutledge D1 (2005) New Zealand's remaining 
indigenous cover: recent changes and biodiversity protection 
needs. Landcare Research, Dunedin, New Zealand. Available at 
http://www.landcareresearch.co.nz/databases/lenz/downloads/ 
New%20Zealand_indigenous_cover.pdf. Accessed 20 September 
2011 


Springer 




296 


Int J Life Cycle Assess (2014) 19:285-296 


Walter C, Stiitzel H (2009) A new method for assessing the sustainabil¬ 
ity of land-use systems (II): evaluating impact indicators. Ecol 
Econ 68:1288-1300 

Weidema BP, Lindeijer L (2001) Physical impacts of land-use in prod¬ 
uct life cycle assessment: final report of the Eurenviro-LCAGAPS 
sub-project on land-use. Department of Manufacturing Engineer¬ 
ing and Management, Lyngby 


WWF, World Wildlife Fund (2012) WWF WildFinder. http://gis. 

wwfus.org/wildfinder/. Accessed 18 August 2012 
ZESPRI (2010) ZESPRI web pages. Zespri kiwifruit Annual Review 
2009/2010, Investors Publications, www.zespri.com. Accessed 6 
July 2012 

ZESPRI (2012) ZESPRI web pages, www.zespri.com. Accessed 6 July 
2012 


4/) Springer 



